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Finite Element Analysis of Axisymmetric Radomes 

Richard K. Gordon, Member, IEEE, and Raj Mittra, Fellow, IEEE 



Abstract — A two-step technique for the analysis (if axisymmet- 
ric radomes is presented in this paper. Initially, an axisymmetric 
Unite element approach is employed, together with an absorbing 
boundary condition for mesh truncation, to determine the near 
fields scattered by an empty radome illuminated by a distant 
source. Next, the reciprocity theorem Is invoked to determine 
the far-field pattern of an antenna encased by the radome, by 
computing the reaction between the current distribution on the 
antenna and the near-field data determined in the first step. The 
details of the formulation are presented along with numerical 
results for two different arrays enclosed by radomes of varying 
permittiyities. 

I. Introduction 

THE problem of determining the modification of the ra- 
diation pattern of an antenna covered by a radome has 
been addressed recently by a number of researchers employing 
a variety of different approaches. Orta et al. [1] have used 
ray optical techniques in conjunction with the reciprocity 
theorem to address the radome problem, and a geometric 
theory of diffraction (GTD) approach has been employed by 
Chikaoka et al. [2] to solve the same. Arvas et al. [3] have 
presented a three-dimensional method of moments solution 
based on the use of the surface equivalence principle, whereas 
Robinson and his colleagues [4 J have used a two-dimensional, 
oblique-incidence, method-of-moments solution to handle the 
problem. A two-dimensional finite element analysis that uses 
an absorbing boundary condition for mesh truncation has been 
employed by Povinelli and D* Angelo [5], and Chan and Mittra 
[6] have developed a spectral-domain approach for the analysis 
of a planar FSS radome. 

In this paper, we present a general -purpose method for ana- 
lyzing axisymmetric radomes with arbitrary inhomogeneities, 
which cannot be conveniently or accurately investigated by the 
ray optical techniques, the method of moments, or approximate 
two-dimensional approaches. We use the reciprocity theorem 
in conjunction with an axisymmetric finite element approach 
that employs a novel absorbing boundary condition for mesh 
truncation. We use a two-step procedure which initially em- 
ploys the finite element method to determine the near fields 
within the radome when it is illuminated by a distant source. 
In the second step, the reciprocity theorem is used to find the 
far -field pattern of an antenna with a given current distribution 
located inside the radome, from the near-field information 
determined earlier in the first step. The method presented here 
assumes that the presence of the radome does not have a 
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significant effect on the antenna currents; it may not work 
if these currents arc significantly changed by the presence of 
the radome. 

II. Formulation 

The details of the two-step procedure for radome analysis 
will now be described. To accomplish the first step, which 
consists of finding the near fields within the radome when it 
is illuminated by a distant source, we employ a finite element 
formulation based on the use of scaled versions of the coupled 
azimuthal potentials (CAP's) 17]. We now present a brief 
review of the approach used in the first step; for further details, 
the reader is referred to [8]. 

We begin by expressing E and H in terms of Fourier series 
expansions in $ as 



fi(p,<M)= £ 



Ez t Tn 



H, 



(1) 



Hz.rn 



(2) 



The potentials we use in this problem are u Tn and v m , where 

u in = E^ m (3) 



(4) 



Inserting (3) and (4) into (1) and (2), and substituting the 

resulting equations into Maxwell's two curl equations, we 

obtain six scalar equations for each value of m. These six 
equations can be written as 



Ep.m = j/m ( m d ^ + k U P fir 



fl.m — 



d{(>v„ 



r dz ) 



jfr, 



( dv m 



+ k 0 pi r 



Op , 



(5) 



(6) 



(7) 



(8) 
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Magnitude of E-theta in the X-Z Plane 



Length of ouiside edge: 2.31 X 



Thickness: 0.1X 
Fig. 2. Radome geometry. 



« m t r + V.(/ m /)f r V(/m m )) + (m/feo)V.(/ ra ^ V(/M; m )) = 0 

(9) 

^n/V+V-(/ m p Mr V(pv m ))-(m/fc o )V-(/ m 0xV(pu m )) = 0, 

(10) 

where 

fm - (fcoMp, z)t r (p, z)p 2 ~ m 2 )~K 

Equations (3H8) demonstrate that £; rrn , E^ m9 E z , m9 
H rt m, Hfam, and #z,m can be expressed in terms of u m 
and This implies that the BOR problem can be solved 
in terms of just the two unknowns, u m and v m , rather than all 
six components of E and H . Furthermore, because of the az- 
imuthal symmetry of the body, rather than a three-dimensional 
mesh, this method requires only a two-dimensional mesh in 
any 4> = constant plane, such as the one shown in Fig. 1. These 
two factors result in a significant reduction in the number of 
unknowns which, in turn, permits the analysis of significantly 
larger radomes. 

Equations (9) and (10) describe the behavior of u m and 
v m in the meshed region. When each of these equations is 
multiplied by a testing function T, and then integrated over 
the entire meshed region T, the resulting equations can be 




•100 -50 o 50 100 

Angle (with respect to the +x axis) 

Fig. 4. Comparison of numerically and analytically determined patterns for 
transparent radome. 

written, with the help of the divergence theorem, as 
J Tu m e r dA- J J f m pt r VT - V(/m m ) dA 

- (m/k 0 ) j j f m VT-(4>x V(/w m )) dA 

= -[ Tf m p, r ^^ldl~(m/k,) 
Jev On • 



Jor 



/J»xV(pt) m )) -ndl 



(11) 



j ^ Tv m/ tr dA ~ J J r frnPVrVT * V(pV m ) dA 

+ (m/ko) J J f m VT • x V(pu m )) dA 

= -/ Tf mP fr%B^<U + (m/ko) 
Jar on 



Jor 



f m T(4>x V(pu m )) ■ fi dl } 



(12) 
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transparent radome. 



are evaluated with the assistance of an absorbing boundary 
condition derived from Wilcox's expansion for the scattered 
fields [9]. This technique is described in detail in [81. 

The finite clement solution for the near fields in the vicinity 
of the radome, when it is illuminated by a plane wave incident 
at some angle 0 with respect to the 2 axis, completes the first 
step of the two-step procedure. Next, the reciprocity theorem 
is used to find the far field at the desired observation angle, for 
any antenna located anywhere within the radome. This is the 
second step of the two-step procedure, and is detailed below. 

Let (■/", M a ) and ( J*, M b ) be two sets of sources existing 



in the same linear medium and let E a and H a be the 
fields produced by (J a ; M a ) alone, and E b and H h be the 
corresponding fields produced by (J b , M b ) alone. Then, by 
invoking the reciprocity theorem, we can show that 



-100 -50 0 50 LOO 

Angle (with retpect to tbe +i axis) 

(0) 

Fig. 5. Effect of radome on antenna pattern for e r = (a) 4.0, (b) 8.0, 
and (c) (4.0. -2.0). 



where dT is the boundary of the meshed region, T. These are 
the equations that are implemented using the finite element 
method. The first term on the right-hand side of each of these 
equations involves an integral of the normal derivative of 
an unknown along the outer boundary dT. These integrals 



(13) 

Thus, for instance, if we set M a and M b to be zero, J a 
to be the electric current distribution on the antenna, and J 6 
to be an electric dipole locaied far from the radome, then 
(13) can be used to determine the far- zone electric field, 
£ a , of the antenna radiating within the radome. In order to 
carry out this calculation, both J a and E h must be known. 
We assume that the presence of die radome has little effect 
on the antenna current distribution; consequently, J a can be 
determined through an analysis of the antenna radiating in 
free spa^e. We also note that E b is the near electric field in 
the vicinity of the radome when it is illuminated by the electric 
dipole J b . But this is precisely what we determined in the first 
step; that is, the plane wave that is used as the incident field 
in the first step is chosen to be the field that would be incident 
upon the radome if it were radiated by J b . Thus, the electric 
field determined in the first step is identical to So, once 
the first step has been concluded for a particular location of 
J b y (13) can be used to find, at this same location, the electric 
field radiated by the antenna within the radome. Similarly, the 
far magnetic field can be determined by setting J b equal to 
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zero and letting M h be a magnetic dipole located far from 
the radome. We note at this point that it is unnecessary to 
recalculate the finite element matrix for each new location 
of J h or M b \ these sources affect only the incident field. 
Thus, the amount of computer time required to solve for 90 
different source locations is only twice that required for a 
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Fig. 8. Horizontal array of dipolcs. 
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Fig. 9. Comparison of nuinerictdly and analytically determined patterns for 
transparent radomc. 

single location. Furthermore, we note that while the radomc 
must have azimuthal symmetry, the antenna may have any 
shape and may be located anywhere within the radome. We 
now consider some numerical examples to demonstrate the 
use of this method. 

III. Numerical Results 

In this section we present some numerical results that 
illustrate the use of the method described above to analyze the 
radome geometry shown in Fig. 2. In each case, we consider 
a radomc that has the shape of a truncated cone which is open 
at the bottom, although the method itself is general enough 
to handle any kind of closure of the radome, as long as the 
azimuthal symmetry is maintained. The length of the outside 
edge of the cone is 2.31A 0 ; its thickness is 0.1A 0 . We will 
present the effects of varying the relative permittivity of the 
radome on the radiation patterns of two antenna arrays. 

The first case we consider is shown in Fig. 3. The antenna is 
an array of three uniformly excited .r-directed di poles located 
on the z axis. The distance between consecutive dipoles is 
0.08o4Aq. We use the two-step procedure described above 
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to determine the far-field value of E$ in the plane 0 = 0° 
for four different radome permittivities: e r — 1.0, 4.0. 8.0, 
and (4.0, -2.0). The results for e r = 1.0 are shown in 
Fig. 4. The reason we consider the case e r = 1.0, i.e., the 
case of a transparent radome, is that for this problem an 
analytically determined pattern can be obtained for comparison 



500 



-500 































——Numerical 1 
— — — Analytical 1 







-100 



100 



-50 0 50 

Angle (witb respect to the +x axis) 

Fig- 1 1 . Comparison of numerically and analytically determined patterns for 
transparent radome. 

with the numerical results. As we can see from Fig. 4, in 
which the dotted line represents the numerical results while 
the solid line depicts the analytically determined pattern, the 
comparison is quite favorable. Note that this does not prove 
the validity of this technique but is a good check. Tn Fig. 5, 
we show the numerical results for c r = (a) 4.0, (b) 8.0, and 
(c) (4.0, -2.0). In each case, the doited line represents the 
radiation pattern in the presence of the radome while the solid 
line depicts the pattern obtained in its absence. We note that, 
in each case, the null at 0°, i.e., in the positive x direction, has 
disappeared. The enhancement of the pattern near -90° can 
be attributed to the fact that the radome is open at the bottom. 
We also note the degradation of the pattern for most positive 
angles, i.e., looking through the radome. 

In Figs. 6 and 7 we consider the same vertical array of 
three dipoles; but while the amplitude of the excitation for 
each dipole is still the same, a phase progression has now been 
introduced. The phase of the excitation on the bottom, middle, 
and top dipole, is 0.0°, -27.86°, and -55.72°, respectively. 
Fig. 6 presents a comparison of the numerical and analytically 
determined results for the transparent radome, while Fig, 7 
shows the effect on the radiation pattern of the presence of 
a radome of relative permittivity (a) 4.0, (b) 8.0, and (c) 
(4.0, -2.0). Like Fig. 4 above, Fig. 6 shows good agreement 
between the numerical results and the analytically determined 
pattern for the transparent radome. Fig. 7 shows the same 
general features that were observed in Fig- 5, but in this case 
the phase progression is present. 

The second case we consider is shown in Fig. 8, The antenna 
is a horizontal array of five uniformly excited a: -directed 
dipoles. The distance between consecutive dipoles is 0.247A 0 . 
We again use the two-step procedure to deterrnine the far- field 
value of Eq in the plane 6 — 0° for four different radome 
permittivities: t. r = 1.0, 4.0, 8.0, and (4.0, -2.0). The results 
for e r = 1.0 are shown in Fig. 9. Just as was the case for the 
vertical array, we see that there is good agreement between the 
numerical results and the analytically determined pattern for 
this horizontal array of uniformly excited dipoles. In Fig. 10 
we see the numerical results for e r = (a) 4.0, (b) 8.0, and (c) 
(4.0, -2.0). We again note that as a result of the introduction 
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Effect of radome on antenna pattern for e, = (a) 4 0 (b) 8 0 
and (c) (4.0, -2.0). 



of the radome, the null at 0° has disappeared, the pattern has 
been enhanced near -90°, and the pattern has been diminished 
at angles near +90°. 

The final results we present are for the same horizontal array 
of dipoles when the amplitude of the excitations is uniform but 
a phase progression is used. Proceeding from left to right, the 



phases of the excitation on the first, second, third, fourth, and 
fifth dipoles are -75.16°, -37.58°, 0.0°, 37.58°, and 75,16°, 
respectively. Fig. J 1 presents a comparison of the numerical 
and analytically determined results for the transparent radome, 
while Fig. 12 shows the effect on the radiation pattern of 
the presence of a radome of relative permittivity (a) 4.0, (b) 
8.0, and (c) (4.0, -2.0). We again note, in Fig. 11, the good 
agreement between the numerical results and the analytically 
determined pattern for the transparent radome. Fig. 12 shows 
the expected effects of the presence of the radome on the 
antenna pattern. 

IV. Conclusions 
In this paper, wc have considered a two-step technique 
for the analysis of inhomogeneous axisymmetric radomes. 
We have verified the approach by computing the radiation 
pattern of an unphased vertical array of dipoles, a phased 
vertical array of dipoles, an unphased horizontal array of 
dipoles, and a phased horizontal array of dipoles within a 
transparent radome for which the radiation pattern can be 
determined analytically. Next, we investigated the pattern of 
these same antennas covered by two lossless radomes and one 
lossy radome. By comparing the results for each case, it is 
possible to determine how the radiation pattern is affected by 
a change in the composition of the radome. Similarly, this 
method could also be used to determine the effect of changing 
the shape of the radome or the addition of a metallic tip. 
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field. The reason for this particular choice of coordinates is 
twofold. First, the use of this normalized coordinate system 
yields very compact equation presentations, and second, these 
coordinates form a two-dimensional Cartesian system in the 
spatial solution domain for the CAP formulation. This latter 
property is of particular value for use with the finite-element 
method, where numerical surface integrations must be per- 
formed over the two-dimensional spatial solution domain. The 
relative constitutive parameters of the medium, e r (R, Z) and 
li r (R, Z), are invariant to the 0-coordinate. 

The equations of the CAP formulation are obtained by first 
decomposing the total fields into azimuthal modes via an ex- 
ponential Fourier series in the 0-coordinate 



E(R t Z t <t>)= ^ e m (R f Z)Qxp(jm<t>) 

m -_oo 
oo 

ri 0 H(R t Z y <t>)= £ *m(K.Z)exp(/m0) 



(la) 



(lb) 



where n 0 = 1 207T £2. Due to the assumed axial symmetry of 
the medium, the modal fields have time-average power orthog- 
onality and decouple when the field expansions in (1) are sub- 
stituted into Maxwell's time-harmonic source-free equations, 
yielding the first-order coupled system given by 



r3<*vm> 



bz 



- jmh 



Z,m 



B(Re, 



0,m J 



a az 



Re r e 



T , 8 «Vm> 1 

T h «- aT"J 

e R,mJ 



(2a) 



(2b) 



(2c) 



(2d) 



(2e) 



(2f) 



where the azimuthal dependence of the modal fields is sup- 
pressed. Using judicious substitution between (2c) and (2f), 
and between (2d) and (2e), it is easily found that all modal 
field components can be generated from two modal scalar po- 
tential functions ^(K, Z, m) and \I/ 2 (R, Z, m), via 



0Xi m =/7 m (*i0X - i?JU r Vt/» 2 ) 



(3a) 
(3b) 
(3c) 
(3d) 



where $ is the unit azimuthal vector, the two-dimensional gra- 
dient operator is defined by 



.9 .3 

V = R ±Z — 

ZR dZ 



and 



f m {R, Z) = ifX r (R t Z)e r (R 9 Z)R 2 - m 2 ]" 1 



(4a) 



(4b) 



As was previously mentioned, the normalized cylindrical 
coordinates (R 7 Z, (p) were chosen only for the practical rea- 
sons of compact presentation and utility in finite-element 
numerical computations. The equations of the CAP formula* 
tion can be easily obtained in an arbitrary coordinate system 
(«1, <M, where u x and u 2 are both orthogonal to 0 and 
form some single-valued system of coordinates in any con- 
stant azimuth planar cross section of the region of revolution. 
For example, the equations throughout this paper can be 
trivially recast into standard spherical coordinates (r, 6, <j>) 
using the simple substitutions R = k 0 r sin 0, Z = k 0 r cos 8 t 
and V = Kl/*o)Wr) + $(l/r)(3/30)]. 

A very important property of the potentials can be ob- 
tained easily by noting that since the modaj^ potentials are 
proportional to the 0-cornponents of e m and /i m , as shown in 
(3b) and (3d), they are continuous everywhere, including at 
dielectric and magnetic interfaces. This property of uniform 
field continuity is very desirable in numerical computations. 
For the case of the variational finite-element method the util- 
ity of the potential field continuity is particularly pronounced, 
in that material interfaces require no supplemental boundary 
conditions and are handled naturally with the same general 
algorithm as is continuously inhomogeneous media. 

III. DIFFERENTIAL EQUATION REPRESENTATION 

The partial differential equations obeyed by the coupled 
azimuthal potentials tf^ and ^ 2 ma Y be obtained directly 
from Maxwell's equations by substitution of the modal field 
generating equations in (3) into (2a) and (2b). The potentials 
are found to satisfy a system of coupled second-order linear 
partial differential equations given by 

V ■ [f m (Re r V\p x + m4> X V^ 2 )] + € r Jp x /R = 0 (5a) 

V ■ l/ m (/*/i,V*2 " ™$ X VlM] + MtzlR = 0 (5b) 

where the gradient operator V is as defined previously Ln (4a) 
and f m (R t Z) is given by (4b). 

To utilize the unimoment technique in obtaining solutions 
to open region radiation and scattering problems requires the 
ability to solve interior Dirichlet type boundary value prob- 
lems within a closed region containing the inhomogeneities of 
the problem, (e.g., an antenna structure or scattering bodies). 
In solving such an interior Dirichlet problem the two-dimen- 
sional solution domain of (5) is a planar constant-^ cross sec- 
tion of an ulterior volume of revolution, as is illustrated in 
Fig. 2. The interior region S, which is bounded by the curve C, 
contains inhomogeneous material specified by the relative con- 
stitutive parameters e r (R* Z) and &(R, Z). The surface S 
could, for example, be the semicircular cross section of a 
geometrical sphere which contains an arbitrarily shaped, (and 
possibly inhomogeneous), material body of revolution, such as 
a torus. 
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Previous chapter 

Jets-In-Crossflow Mixing 

Clifford E. Smith, Principal Investigator 
Co-investigators: Daniel B. Bain and James D. Holdeman 
CFD Research Corporation/NASA Lewis Research Center 

Research Objective 

The objectives of this work were to identify improved mixing schemes for rich-burn/quick-mix/lean- 
burn combustors applicable to advanced aircraft engines and to assess the effects of design parameters 
on mixing effectiveness. Efforts were focused on jet mixing in rectangular cross-sectional geometries. 

Approach 

Numerical parametric studies were performed on different orifice configurations. The computational 
analyses were performed using a three-dimensional Navier-Stokes flow solver (CFD-ACE) capable of 
analyzing turbulent reacting flows in complex geometries. The numerical results were examined using 
an interactive visualization package (CFD-VIEW). 

Accomplishment Description 

Over 60 orifice configurations were analyzed and assessed. The modeled geometry consisted of a 
rectangular duct with a row of jets on the top and bottom walls. Isothermal flow conditions were 
assumed. Systematic, parametric analyses were performed, consisting of variations in (1) orifice lateral 
alignment (in-line and staggered), (2) jet-to-mainstream (J) momentum flux ratio (16, 36, 64), (3) orifice 
aspect ratio (4:1, 2:1, 1:1, circle), (4) jet-to-mainstream mass flow (MR) ratio (2:0, 0.50, 0.25), and (5) 
orifice spacing-to-duct height (S/H) ratio (0.125 <= S/H <= 1.5). The numerical calculations were 
performed with grids varying in size from 50,000-150,000 cells that required approximately 2-4 Cray Y- 
MP hours and 10-15 megawords of memory. 

Significance 

At optimum S/H, in-line lateral arrangements produced faster initial mixing than non-impinging, 
staggered arrangements due to their smaller geometric orifice size. Previous publications show that S/H 
((root) J) is a practical design parameter for mixing optimization even at high MR, and increasing J 
improves initial mixing at optimum S/H. 

Future Plans 

The effect of mass flow ratio and aspect ratio on mixing effectiveness will be assessed, and emission 
characteristics will be inferred from existing cold-flow data. 
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Publication 

Bain, D. B.; Smith, C. E.; and Holdeman, J. D.: CFD Mixing Analysis of Axially Opposed Rows of Jets 
Injected Into Confined Crossflow. AIAA Paper 93-2044, June 1993. 




Typical jet mixing numerical results. In-line circular orifices: S/H = 0.375, J = 36, MR = 2.0. Jet mass 
fraction: red = 0.0 to blue = 1 .0. 

Next chapter 
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